o 
o 

(N 



S3 

c 

I 

C/3 



I 

C 

o 
o 



> 

in 
m 
in 
o 

o 

c3 



i 

c 

o 
o 



X 
S3 



New Real space method for calculation of physical properties of a disordered system. 
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We introduce a new real space super cell approximation method for treating the electronic states of 
disordered systems. This method is general and allows both randomness in the on-site energies and 
in the hopping integrals. In the special case of randomness in the on-site energies only, this method 
is equivalent to the Non Local Coherent Potential Approximation (NLCPA) derived previously. 



For the past three decades the Coherent Potential 
Approximation (CPA) was the most widely used self- 
consistent mean field treatment of disordered systems 
As a single site approximation, the CPA Green's func- 
tions are averaged over the randomness on a single site, 
and so could not handle the effects of correlations be- 
tween the random potentials on neighbouring sites. In 
order to include such inter site effects extensions of 
CPA were developed using the locator formalism 0, 
However, recently the Dynamical Cluster Approximation 
(DC A) has been introduced for interacting systems 0. A 
version of DCA has also been used to treat disordered 
systems [(§■ This method extends CPA by including a 
weakly wave vector dependent self energy E(fc;w). The 
DCA was originally derived on basis of restricting the 
role of momentum conservation in calculating the self- 
energy, by applying a coarse graining procedure to the 
First Brillouin Zone (FBZ). Momentum conservation is 
applied only to the coarse grained wave vector K„ in the 
FBZ and in the corresponding Laue function. In this 
method a key computational step makes use of a Fourier 
like transformation of coarse grained wave vectors K„ to 
a real-space set of coordinates, R„. These coordinates 
were interpreted as the lattice sites of a real cluster of a 
real lattice, however this connection was not established 
directly. 

By using effective medium theory in real-space we de- 
veloped a non-local CPA (NLCPA) method. In this 
method a cluster of impurities is embedded in an effective 
medium, in such a way that all cluster sites are equiv- 
alent and lattice periodicity is maintained for the dis- 
order averaged Green functions. For a one dimensional 
binary alloy numerically we showed that a periodic super 
cell density of states is similar to the NLCPA density of 
states|3- We also applied NLCPA as a real space super 
cell approximation to explain the resonance peak that 
appears in the density of states of substituted Zn im- 
purity in CuOi plane of BiiSr-iCaiCu^O^i,^- Subse- 
quently NLCPA was also extended to the case of the first- 
principles Korringa-Kohn-Rostoker (KKR) CPA model 
for the electronic structure of disordered systems. 

In this paper we show that the super cell self energy 
periodicity in real space leads naturally to coarse graining 
in k-space, providing us with an alternative and general 



derivation of the NLCPA method. This establishes di- 
rectly the connection between the real-space cluster and 
the coarse graining in k-space, showing clearly for the 
first time that both DCA and NLCPA are equivalent to 
a super cell approximation in real space. 

We start our investigation on a general tight binding 
model for a non-interacting alloy system, 



(1) 



where ty are the hopping integrals which may be random, 
fj, is the chemical potential and £j is the random on site 
energy which takes values —5/2 with probability 1 — c for 
host sites, and 5/2 with probability c for impurity sites. 

The equation of motion for electrons with the above 
Hamiltonian, Eq^ for any general impurity configura- 
tion {e{\ in terms of the Green function is, 

( (w - + + tu ) G{1, j;u) = S i:j (2) 

i 

We expand the full random Green function G(i,j;uj) 
in Eq|2 in terms of the clean system Green function 

G(i, j;u>) = G°(i, + £ G°(i, I; u)V u > G{1 (3) 

«' 



where the random potential matrix V u > is 
v w = £ i s u' + §t u' > 



(4) 



5t {l > — t u i —t% is the difference between hopping integral 
in random and clean system and G (i,j;u>) is 

G°(i, j; u) = ^J2 e ^ 3 {u-^ + nY 1 . (5) 



where e k = -^Y,, 1 ,,■ 



„k.i 



is the band structure for the 



clean system. 

The Dyson equation corresponding to Eq[3]is 

G(i,j;u) = G (i,j;oj) + Y / G°(i,1;lo)1:(IJ';lo)G(1' ,j» 
U' 

(6) 



2 



where, 

J2(V ll >G(l\j;w)} = J2X(lj'^)G(l',j;u). (7) 
l' i' 

The disorder average Green function, EqEl can be writ- 
ten, 

G(i, = 1 £ ^ 3 ( ^ - e k + M - E(k; w) ) _1 . 

k 

(8) 

where 

S(k;o;) = l^e- sk - r «E(i I i;a;). (9) 

is the Fourier transform of the self-energy. 

Fourier transformations of this system of equations 
obey the orthogonality condition, 

k 

where indices i and j run over all N lattice sites in the 
crystal. Because the exact solution of Eqs|3| and [H is 
impossible, they are usually solved in different single site 
approximations such as the Born approximation, the T- 
matrix approximation and the CPA.|7j It is well known 
that in all these approximations inter site impurity effects 
has been neglected. 

Here we introduce a new approximation, and we show 
that, for the case of St H ' = 0, it is equivalent to the 
NLCPA §. The derivation of the approximation is as 
follows. Consider a random alloy system that has been 
divided to iV c super cells. The super cell Green function 
G SC (I, J; u>) can be expand as 

G ac (I,J;u) = G°(J,J;w) 

+ 5]G (/ ) iiw)VB'G,c(i' ) J;w). 

i,v 

(11) 

where I and J refer to sites in the same super cell and 
1,1 indices refer to the whole lattice. The average Green 
function corresponding to Ea llll is. 

G sc (I,J;w) = G°{I,J;w) 

+ ]T G°(I, I; w)£ sc (Z, I'; u>)G sc {l', J; w). 
l,l' 

(12) 

Now we apply an approximation that is based on 
two assumptions: first neglect the hopping integrals 
deviation, 8t u > , between impurity sites of different super 
cells, that is, 8t u i = when I and I are in different su- 
per cells and second neglect multiple scattering between 
impurities in different super cells. Therefore we have 

E sc (i,j; u>) = 0, if i and j (£ same super cell. (13) 




'(a) ' '(b) ' 



FIG. 1: (a): A binary alloy divided into super cells containing 
four sites, (b) The approximation of impurity averaging over 
the super cells. By neglecting the multiple scattering between 
impurities in different super cells, the self energies in different 
super cells obey the periodicity law given in Eq|14l 
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FIG. 2: A two dimensional lattice divided into super cells 
of nine sites. The arrow indicates the super cell periodicity 
vector rjv e = 3ai + 3a2, where Nd = 3 and N C 2 = 3. 

Figffl shows a binary alloy system divided into periodic 
super cells containing four impurity sites. By neglecting 
both the hopping integrals deviations 8t U ' and also the 
correlations between different super cells sites, the self en- 
ergy in each super cell is independent of the other super 
cells and is periodic with respect to super cell transla- 
tion vectors. These conditions imply that the disorder 
averaged system self energy has the following super cell 
periodicity, 

S sc (r/j + r Nc ; ui) = S sc (r/j; u>) (14) 

where I and J refer to the sites in the same super cell, 
and the super cell translation vector tn c is 

3 

rjv c = ^2 mi N ci a l; (15) 

i=l 

where N c — iV cl N c2 N c3 is the number of lattice sites 
in a 3-dimensional super cell, a 4 are primitive vectors of 
real lattice, and N ci is the number of sites in the super 
cell in the ith direction and mj are integers. Fig|2| shows 
r^v c for a two dimensional lattice with periodicity over 
nine sites in each super cell. Note that when iV c — * oo 



3 



all super cells will coincide and are equal to the full real 
lattice. 

Fourier transformation of Eg 1 141 implies that we should 
restrict the self energy wave vectors k to those obeying 



= ik.rjv c _ 



1. 



These are the set 



(16) 



(17) 



: i 
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; i 





where bi are the reciprocal lattice primitive vectors and 
li are integers. Notice that number of the K„ in the first 
Brillouin zone is equal to number of the lattice sites, N c , 
in a super cell. Therefore the Fourier transformation of 
the self energy is 

S(K„; W ) =l^e- lK - r y£ sc (i,j;u;). (18) 



1,3 



By inserting Eq^Hinto Eq^Jwe find that 

E(K n ; u) = ±J2 ^ lK - riJ Zsc(I, J; to) (19) 



i, J 



where / and J refer to the sites in the same super cell. 
By converting the Eq^|to real space we find 

S sc (7, J; u) = ±r e- K - r "S(K n ; w) (20) 



K, 



where the orthogonality condition in the super-cell is, 

i_5>-*K„. r „ =( ^ (21) 



K„ 



Because the set of {K„} vectors divides the volume of 
the first Brillouin zone to Nc equal patches, we therefore 
identify each patch, n, by the corresponding wave vector 
K n that is located at its center. Therefore we define the 
relation between K„ and k inside each patch as follows 



k = K„ + k„ 



(22) 



where k n are the wave vectors inside of nth patch with 
respect to the center of the patch. Fig|3] illustrates this 
relationship, Eq[53 for one of the patches in a two di- 
mensional system for N c — 9. By comparing Eas. l21l and 
HOI we find that the super cell approximation is equivalent 
to the replacement 



(23) 



Considering of the Eqs. EU Eld and Eq.[Hwe nn d that 
the super cell approximation Green function G sp (I, J; lo) 
is 

G ac (J,J;w) = -^^e* K - r «G(K n ;u/) (24) 



FIG. 3: Relation between k, K n and k' n in the First Brillouin 
Zone (FBZ) for a two dimensional system with a 9 site super 
cell. As given by Eqs. 1171 and 1221 there are 9 different K n in 
the FBZ, and so the FBZ is divided in to 9 equal areas, while 
each K„ is located at the center of one of these areas. 



where 



N 



X] ( u ~ e K„+k'„ + M - S(K„; u) 

k' 



(25) 

In short Eq. ^| and 1231 are the super cell approximation 
conditions. In the limit of N c — * oo, k^ — > and also 
K„ — > k. Therefore Eas. 1211 and 1191 respectively convert 
to Eas. 1101 andl9lin the exact average system. Also for the 
exact average system in the limits of super cell periodicity 
with N = mN Cl with m an integer, in the case where 
correlation between super cells is neglected, the exact 
average self energy with these conditions equals the super 
cell periodic self energy, that is 



S(J, J; a) = E ac (J, J;w). 



(26) 



By use of Eq. 123 the Fourier transformation of Eq. 1261 
can be carried out by reordering the summation over k 
as J2k = Dk„ J2k' n and usin g the fact that 1 = ik- 
We find that 



S(K n ;a;) = ^^S(K n +k^w) 

C k' 



(27) 



and Eqs. E3 El and El will become equal to Eas. I2"T1 [271 
and 1241 respectively. This means that the super cell cal- 
culation in an effective medium is an approximation for 
the exact disordered system that includes only the cor- 
relation between impurities within a single super cell. 

On the other hand, by eliminating the clean Green 
function matrix Go between Eqs. |21 and the random 
matrix Green function G can be expressed in terms of 
the average Green function G as 



G = G + G(V - S)G. 



(28) 



K„ 



By applying the super cell approximation with Eq. 1131 
and taking the average over impurities in all super cells 
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FIG. 4: A four site impurity super cell in an effective medium 
of mean-field super cells. 

except one, Eq. |^H1 reduces to the following equation for 
the remaining super cell, 

G™^ - G sc + G SC (V SC - S sc )G™p. (29) 

where all the matrices in Eq.[2H|are of dimension N c x N c . 
By rearranging Ea l29l we have 

G^" 1 +V SC = G- 1 + S SC . (30) 

Defining each side of Eq[3U| as a cavity super cell Green 
function Q we have 

G imp-1 +Vsc =g- 1 . (31) 

and 

G" 1 + £ sc = (32) 

Fig. 0] shows an example four site super cell approxima- 
tion in which the average has been taken over all super 
cells except one impurity super cell. The matrix elements 
of Eq.EJare 

Gl™P(I,J-u;) = Q{I,J;w) 

+ £$(I,L; W )y^G™*(i/,.7;w). 

L,L' 

(33) 

where the cavity Green function £/(/, J;u>) obeys 

Q(I,J;u) = ^- c J2^- riJ G(K n ;u;). (34) 

Similarly, the matrix elements corresponding to Eq. 1321 
are, 

G. c (J,J;w) = Q{I,J;u) 

+ J2 L 5 w)S sc (L, L'; oj)G sc {L', J; u). 

(35) 

where the average of the impurity Green function 
G l ™ p (I, J;uj) over all impurity configurations {V/j} in 
the super cell is 

(G i ™v(I,J;u)) = G sc (I,J;u). (36) 



Eqs. I2SJ |23 EUl EH constitute a closed set to be solved 
self consistently. 

The above system of equations can be implemented 
numerically by the following algorithm: 

1- make a guess for E(K„;w), usually zero. 

2- Calculate G(K n ;w) from Eq.^SJ 

3- Use Eq. to calculate the Fourier transform of the 
cavity Green function 

g(K n ;uj) = (G- 1 (K n ;o;) + E sp (K n ;a;))- 1 (37) 

4- Calculate the impurity super cell Green function 
G™ p (/, J;uj) from Eq.ESJ 

5- Calculate the average Green function from Eq. 1361 and 
Fourier transform it to K„ space by using Eq. |^ 

6- Using Eq.|2Z|calculate the new self energies E SC (K„; u) 
and go back to step 1 and repeat whole process until con- 
vergence has been obtained to a desired accuracy. 

As an application of this method we calculate the den- 
sity of states of a two dimensional square lattice at half 
band filling, n = 1 and S = 6t° in which N c = 1 (CPA) 
and N c = 9. For the N c = 9 case, we also show the effects 
of introducing random hopping parameters, by consider- 
ing the two cases, StAA = 0, StAB = and StAA — 4i°, 
StAB = t , where 

— *AA ~ t (38) 

which i and j are nearest neighbour sites and both are A 
type, and we also defined, 

St {ij) = t AB - t° (39) 

where at the nearest neighbours sites i and j one is A 
and another is B type atom. Fig|5] shows that for the 
Nc = 1 (CPA) and the case of Nc = 9 without random 
hopping, in which S = 6t°, we are at the band splitting 
regime [ill |l2|. In this case at half band filling a metal- 
insulator phase transition is taken place, in spite of dif- 
ferent gap sizes. While for the case of Nc = 9 with ran- 
dom hopping StAA = 4t°, StAB — t°, the band splitting- 
is not happened, thus the system is a metal. Therefore 
due to including randomness in both the on-site energies 
and also the hopping integrals, EMSCA technique can be 
provide a more realistic results. 

Finally, we summarize our discussion with the follow- 
ing conclusions. We have introduced a new real space 
effective medium super cell approximation (EMSCA) for 
random systems. For the case N c = 1 this recovers the 
single site CPA formalism. However for larger values of 
N c the method extends CPA by including both the ef- 
fects of multiple scattering, and by allowing randomness 
in the hopping integrals. Furthermore, in the limit of 
N c — > oo the method is exact. The EMSCA in the spe- 
cial case where hopping randomness is neglected, St = 0, 
leads to an alternative derivation of the NLCPA tech- 
nique for disordered systems. This derivation completely 
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FIG. 5: The density of states of a two dimensional square 
lattice of a binary alloy A c Bi- c at half band filling, n = 1, 
for c = 0.5 and 8 = 6t° when iV c = I (CPA) and N c = g . We 
compared N c — 1 DOS and of two cases of N c = 9, first 
StAA = 0, StAB = and second StAA = it , StAB = t°. For 
both non hopping integral randomness cases the band split- 
ting happened ( metal insulator phase transition take placed), 
while for the other not (still is a metal). 



establishes the NLCPA as a valid and useful extension 
of the old and popular CPA method, which incorporates 
the effects of inter-site correlations. 
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